MZ-TH/09-18 
larXiv: 0905 .4 465 [hep-ph] 



0(a s ) corrections to the polar angle 
dependence of the longitudinal spin-spin 
correlation asymmetry in e + e~ —> qq 

S. Groote, J.G. Korner and J. A. Leyva 

Institut fur Physik, Johannes-Gutenberg-Universitat, 
Staudinger Weg 7, 55099 Mainz, Germany 



Abstract 

We provide analytical results for the 0(a s ) corrections to the polar angle dependence 
of the longitudinal spin-spin correlation asymmetry in e + e~ — * qq. For top quark 
pair production the 0(a s ) corrections to the longitudinal spin-spin asymmetry are 
strongly polar angle dependent and can amount up to ~ 4% in the t/ 2 -range from 
above ti-threshold up to \[q* = 1000 GeV. The 0(a s ) radiative corrections to the 
correlation asymmetry are below ~ 1% in the forward direction where the cross 
section is largest. In the e + e~ — > bb case the 0(a s ) corrections reduce the asymmetry 
value from its = value of —100% to approximately —96% for g 2 -values around 
the Z peak and are practically independent of the value of the polar angle 9. This 
reduction can be traced to finite anomalous contributions from residual mass effects 
which survive the mj, — ► limit. We discuss the role of the anomalous contributions 
and the pattern of how they contribute to spin- flip and non-flip terms. 



1 Introduction 



Since polarized top and antitop quarks decay weakly before hadronizing they retain their 
original polarization from the production process when they decay. There are thus signifi- 
cant correlations between the decay products of jointly produced top and antitop quarks. 
The relevant spin information is encoded in the spin-spin density matrix of the (tf)-system. 
For hadronically produced (tt)-pairs such correlations have been thoroughly discussed in 
the review article [lj. In e + e~ — > ti interactions the spin-spin density matrix of the (tt)- 
system (2J [3J HI El EJ [7] and the ensuing correlations between the decay products [21 [3] have 
been studied in a number of papers. 

Of interest is also the role of quark mass effects in the production of quarks and gluons in 
e + e~-annihilations. Jet definition schemes, event shape variables, heavy flavour momentum 
correlations [HI 0, EE] and the polarization of the gluon [H] are affected by the presence of 
quark masses for charm and bottom quarks even when they are produced at the scale of 
the Z° mass. A careful investigation of quark mass effects in e + e~-annihilations may even 
lead to an alternative determination of the quark mass values [8j [91 [101 [12] . There is an 
obvious interest in quark mass effects for (tt)-production where quark mass effects cannot 
be neglected in the envisaged range of energies to be covered by the proposed International 
Linear Collider (ILC). At the low mass end quark mass effects are important in the m q — > 
calculation of radiative corrections to quark polarization variables because residual mass 
effects change the naive non-flip pattern of the m q = polarization results [131 EH]- I n 
QED the same 0(a/ir) residual mass effects have been widely discussed in polarized lepton 
production [T5l [16], [T7] and polarization effects in lepton decays [HI [191 120] - 

In this report we provide analytical results for the 0(a s ) radiative corrections to 
the longitudinal spin-spin correlation asymmetry (called longitudinal spin-spin asymme- 
try for short) and its polar angle dependence for massive quark pairs produced in e + e~- 
annihilations. This goes beyond the calculation of Ref. [51 [6] where the polar angle depen- 
dence was averaged over. In contrast to the analysis [7J, which was numerical, our results 
are presented in closed analytic form which allows one to explicitly study the high energy 
(or m q — > 0) limit of the relevant spin-spin density matrix elements which differ from the 
naive m q = values due to the residual quark mass effects. 

The longitudinal polarization of massive quark pairs produced in e + e~-interactions 
affects the shape of the energy spectrum of their secondary decay leptons. For example, 
the longitudinal spin-spin correlation effects in pair produced top quarks and antitop quarks 
will lead to correlation effects of the energy spectra of their secondary decay leptons and 
antileptons. We discuss in detail the m q —>■ limit of the relevant polarized structure 
functions and the role of the 0(a s ) residual mass effects which are mostly relevant for 
charm and bottom quarks for energies at and above the Z and for top quarks at the 
highest energy y/q* = 1000 GeV. We delineate how residual mass effects contribute to 
the various spin-flip and non-flip terms in the m q — > limit for each of the three spin-spin 
correlation functions that describe the polar angle dependence of the longitudinal spin-spin 
correlations. 
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2 Joint quark- ant iquark density matrix 
and polar angle distribution 

Let us begin by defining the differential joint quark- ant iquark density matrix d<sr a = 
^ cr AiA 2 -A'A' where Ai (X[) and A2 (A 2 ) denote the helicities of the quark and antiquark, 
respectively. The label a specifies the polarization of the initial 7* and Z, or interference 
contributions thereof, which determines the polar angle dependence of the cross section and 
the longitudinal spin-spin correlations. In our notation the three polarization components 
are called a = U (unpolarized transverse), L (longitudinal) and F (forward-backward). 

In this paper we concentrate on the longitudinal polarization of the quark and antiquark, 
and in particular, on their longitudinal spin-spin correlations. Thus we specify to the 
diagonal case Ai = \[ and A 2 = A' 2 . The diagonal part of the differential joint density 
matrix can be represented in terms of its components along the products of the unit matrix 
II and the z-components of the Pauli matrix gr 3 (<p~ 3 = p±e for the quark and «r 3 = p 2 e for 
the antiquark, pi = Pi/\f>i\). One has 



where the Pauli matrices to the left and right of the <8> symbol are associated with the 
quark and the antiquark, respectively. 

In order to get a feeling about the physical significance of the single and double spin 
functions at x ^ and o^ 1 ^ consider the case where the polarization of the quarks and 
antiquarks is purely longitudinal in the helicity system. This would be the case in the high 
energy limit since the transverse spin components of the top and antitop quark go to zero 
as m q / \fq*. The differential distribution with respect to the polar angles cos 9i 2 specifying 
the angles between the (quark, antiquark) and given decay products thereof reads 



where the asymmetry parameters ai j2 specify the analyzing power of the particular final 
state in the respective decay channels. Consider, for example, (tf)-production followed by 
the semileptonic decays t — ► b + £ + + vg and t — > b + 1~ + vg. Each of the decay products of 
the top and antitop can serve as a polarization analyzer of the spin of the top and antitop 
quarks. From the experimental point of view the charged leptons £ ± are ideally suited for 
this purpose because they are easy to detect and because their analyzing power is maximal 
with ai j2 = ±1 (see e.g. [2T] ) . 

An alternative but equivalent representation of the longitudinal spin contributions in 
Eq. (00) can be written down for the double density matrix elements a a (s{, s|) in terms of the 
longitudinal spin components s[ = 2Ai and s e 2 = 2X 2 with s[, s 2 = ±1 (or s[, s 2 G {T, ID- 
One has 




(1) 




(2) 




(3) 
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Eq. is easily inverted. The result is 



°« = CTa(TT) + 0"q(T1) + ^a(lt) + ^(11) 

V { a l] = ^(TT)+Cr a (U)-^at)-Cr a (U) 

^ = ^(tT)-^(u)+^at)-cr a (u) 

ft**) = ffa (tt) - a a (U) - Cr a (|t) + Cr a (||) 



(4) 



In the following we shall refer to the (fl) and (||) spin configurations as aligned even if 
the quark and antiquark are not parallel. 

0(a s ) radiative corrections to the unpolarized rate components <y a have been discussed 
before (see e.g. [13j [HJ [22]) including beam polarization effects [H] and beam-event corre- 
lation effects [HI [22]. The 0(a s ) radiative corrections to the longitudinal spin component 
&u+L ( an d thereby cr^ L ) have been calculated in Refs. [HI [23]. Beam polarization and 
0(a s ) beam-event correlation effects for (a = U, L, F) were calculated in [HI [241 [25] . 
As concerns the spin-spin correlation asymmetry a^ 2 ' the 0(a s ) tree graph contribu- 
tions have been written down in Ref. [I]. The complete 0(a s ) radiative corrections to 
the fully integrated spin-spin correlation component crjj^i were calculated in Refs. [5], [6] 
where beam-event correlation effects were averaged over. These results were confirmed in 
a numerical calculation by Brandenburg et al. [7] . As mentioned in the introduction to this 
paper we are aiming to determine the polar angle dependence of the longitudinal spin-spin 
correlation asymmetry, i.e. we are interested in the polar angle structure induced by the 
rate functions a a and the longitudinal spin-spin functions a^ 2 ' (a = U, L, F). 

As before we write the electroweak cross section and the spin-spin correlation compo- 
nents in modular form in terms of three building blocks |14| . namely the lepton tensor 
(which encapsules the beam polarization parameters and the angular dependences), the 
hadron tensor (which contains the hadron dynamics) and an electroweak coupling matrix 
gij(q 2 ) which connects the two. Since the electroweak model dependence and the polar 
angle dependence of a a and da^ 1 ^ are the same, we introduce a compact notation a^ ili2 ^ 
where a^ ili2 ^ stand for either a a or a^ li2 \ Thus we write 

1^9 = s (1 + cos e) \w +9l2au 

3 „ / 3{44} , „ 4{^£ 2 } 



+^cos8{ K g i3 a^+guar i2i )- (5) 

The index j = 1,2,3,4 in a^ lt2 ^ runs over the four linear combinations of bilinear products 
of vector and axial vector currents defined by 

~1 („VV 1 „AA\ 2 ^-(„VV „AA\ 

°a = ^ a + a " a = 2^ a ~ a 

^3 1 („VA „AV\ 4 ^fJVA , „AV\ ( a \ 
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6 is the polar angle between the electron beam direction and the top quark direction. 

The matrix (i,j = 1, 2, 3, 4) specifies the electroweak model dependence of the e + e~ 
cross section. For the present discussion we only need the components gn, gu, 943 and g^. 
They are given by 

9u = Q 2 f -2Q f v e v f Rexz + (vj + a 2 e )(v 2 f + a 2 f )\xz\ 2 , 

9n = Q 2 -2Q f v e v f Rexz + (vj + a 2 e )(v 2 - a 2 )\xz\ 2 , 

943 = 2Q f a e a f lmxz, 

#44 = -2Q f a e a f Rexz + 4:V e a e v f a f \xz\ 2 (7) 

where for the Standard Model Xz{q 2 ) = gM 2 z q 2 /(q 2 - Mf + iM z T z ) with M z and T z 
the mass and width of the Z° and g = G F (8V2na)- 1 m 4.49 x 10- 5 GeV~ 2 . Q f are 
the charges of the final state quarks to which the electroweak currents directly couple, v e 
and a e , Vf and a,f are the electroweak vector and axial vector coupling constants. For 
example, in the Weinberg-Salam model one has v e = —1 + 4 sin 2 9w, &e = — 1 for leptons, 
Vf = 1 — | sin 2 Q w , a>f = 1 for up-type quarks (Qf = |), and Vf = — 1 + | sin 2 W , a,f = — 1 
for down-type quarks (Qf = — |). At low energies, where only QED interactions survive, 
the only remaining components of the electroweak coupling matrix g^ are gn = Q 2 and 
9i2 = Q 2 - In this paper we use Standard Model couplings with sin 2 6 W = 0.231 and 
M z = 91.188 GeV, T z = 2.487 GeV [26]. 

The hadronic building block is entirely determined by the hadron dynamics, i.e. by the 
current-induced production of a quark-antiquark pair which, in the 0(a s ) case, is followed 
by gluon emission. At 0(a s ) one also has to add the one- loop contribution. We shall work 
in terms of unpolarized and polarized hadron tensor components H° a and H^ 1 ^ 2 ' where 
the spin decomposition is defined in complete analogy to Eq. Q. The notation closely 
follows the notation used of Ref. |14j . 

In the two-body case e + e~ — > q(pi)q(p2) (Born and one-loop contribution) the unpo- 
larized rate components a J a and the longitudinal spin-spin components related 
to the corresponding two-body hadronic tensor components H^ il&2 ^ via 

J{W2} = !^m{«2} (o\ 

where v 2 = 1 — Am 2 /q 2 :— 1 — £. The helicity structure functions H^ 1 ^ (a = U,L,F) 
are obtained from the hadron tensor components H^ 1 ^ via covariant projection 

TT3{ht2) _ ( _~ _ PlPl\ TTj{£l£ 2 } TTj{tli2} _ PlPlTT{il£ 2 } 

n U — 9iiv 2 n \w ' L — 9 [a> ) 

\ Plz ) Plz 

H F ^ = i£ ^2l?%Hti M ^ (9) 
PizVQ 

where g^ v = g^ v — q^qv/q 2 and pi^ = p^ t — (pi ■ q)q^/q 2 are the four-transverse metric 
tensor and the four-transverse quark momentum, respectively. The two-body hadron tensor 
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Hj£ li2 ^, finally, is denned by the product of matrix elements 

HiM^s[,s l 2 ) = (g(4)?(4)l^|o)(o|Jr >A k(4)g(4)}, (10) 

where, again, the index j = 1,2,3,4 specifies the current composition in terms of the 
parity-even (for j = 1,2) and parity-odd (j = 3,4) products of the vector and the axial 
vector currents as defined in Eq. (jBJ). 

In the three-body case e + e~ — > q(pi)q(p2)g(P3) the hadron tensor W is obtained from 
the square of the current matrix elements according to 

Hf iu (p u p 2 ,p 3j s[,s l 2 ) = £ (?(4)^l-^ A |0><0|Jj v;A |?( a i)?(4)P>- (11) 

gluon spin 

Note that the three-body hadron tensor in Eq. ( TTTT) has a mass dimension differing from 
that of the two-body hadron tensor in Eq. (TlOT) . The projection of the three-body hadron 
tensor onto the six helicity structure functions H^ il£ ^ (a = U, L, F) is done in identical 
fashion to the two-body case. What is finally needed is the relation between the differential 
rate components and the three-body helicity structure functions. The relation is given by 



1,1 Hi^(y,z)\ (12) 



daiVUV _ ira 2 v 
dydz 3q 4 \ 16n 2 v 

where we have explicitly referred to the phase space dependence of the three-body hadron 
tensor in order to set it aside from the two-body hadron tensor in Eq. ( flOl) . As kinematic 
variables we use the two energy-type variables y = 1 — 2piq/q 2 and z = 1 — 2p 2 q/q 2 . 

We mention that the above formalism allows for a straightforward incorporation of 
transverse and longitudinal beam polarization effects [TJ] . For example, if one starts with 
longitudinally polarized beams, one has to effect the replacement 

g u -> (1 - h-h + )g u + (h- - h + )g 4i (i = 1, 2) 

94j - (l-h-h+)g 4j + (h- -h + )g Xj (j=3,4) (13) 

where h~ and h + (—1 < < +1) denote the longitudinal polarization of the electron and 
the positron beam, respectively. Clearly there is no interaction between the beams when 
h + = h~ = ±1. The additional electroweak components g 4 i, g 42 , gi3, and gi 4 needed in 
Eqs. ( Tl3|) are given by 

041 = 2Q f a e v f Rexz - 2v e a e (v 2 + a 2 )\xz\ 2 , 

042 = 2Q f a e v f Rexz ~ 2v e a e (v 2 - a 2 )\xz\ 2 , 
g 13 = -2Q f v e a f lmxz, 

0u = 2Q f v e a f Re xz ~ 2(i% + a 2 e )v f a f \xz\ 2 ■ (14) 
The incorporation of transverse beam polarization is described in Ref. 
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3 Born term contributions 



Let us begin by listing the nonvanishing Born term contributions to the unpolarized and the 
polarized spin-spin two-body hadron tensor components. They can be obtained from the 
two-body hadron tensors in covariant fashion by using the projections ([9]) and the two-body 
representations of the spin vectors (|25j) and (|26j) . An alternative and rather convenient 
approach is to use helicity amplitudes h^l x A? in the cm. frame where A denotes the rest 
frame m-quantum number of the (j*,Z) (the quark momentum defines the ^-direction). 
The helicity amplitudes are given by 



' 2 2 v ' 2 2 

h v ±l . M , = -V2V^, h*,= T V2v^. (15) 



One obtains 



v 2 ) 



and 



Hlj(Born) = 2N c q 2 (l + v 2 ), 

H 2 (Born) = 2N c q 2 (l-v 2 ), 

Hl(Born) = H 2 (Born) = N c q 2 (l 

H^(Born) = AN c q 2 v (16) 

H^\Born) = -2N c q 2 (l + v 2 ), 

H 2(eie2) {Born) = -2N c q 2 (l - v 2 ), 

Hi} tlta) (Born) = H 2 L {h " 2) (Born) = N c q 2 (l-v 2 ), 

H 4 / 1&2 \Born) = -AN c q 2 v. (17) 

For completeness we also list the single spin Born term contributions. They read 

H 3 /^ = 0, 

= ±AN c q\ 

H 1 /^ = ±2N c q 2 (l + v 2 ), 

H 2 F ih ' 2) = ±2N c q 2 (l-v 2 ), (18) 

where the upper and lower signs stand for the quark and antiquark case, respectively!]] 
Note that one has the Born term relations 

H]j 2 (Born) = -H^^ (Born), 
H^(Born) = -H*/ 1 ^ (Born), 

Hl' 2 (Born) = H 1 /^ (Born). (19) 



1 The signs of H 1 ^ 212 differ from those in Ref. [27] because, contrary to [27], in the present paper we are 
always referring to the same polar angle in Eq. ([5]) for both the quark and antiquark case. 
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These relations can be seen to arise from angular momentum conservation at the Born term 
level where the produced quark and antiquark are in a back-to-back configuration. The 
easiest way to see this is by noting that at the Born level one has Hu,f('\ j \) = Hu,f(11) = 
and i?i(T!) = Hl{[]) = from angular momentum conservation. The above Born term 
relations then follow from inserting the above zero entries into Eqs. (jlj) (with a a replaced 
by H a ). It is quite clear that the relations Eq. ffl~9l) no longer hold true at 0(a s ) in general 
since the quark and antiquark are no longer back-to-back due to additional gluon emission. 

Let us now introduce the normalized cos ^-dependent longitudinal spin-spin asymmetry 
P u (cos9) which is defined by 

P « (cos g] _ l(l + c OB ^) g r" + jBin^ g ^ + |co s 9^"» 

1 ; ~ §(l + co S 29)<7 [/ + fsin 2 9<T i + fcos0<7 r 1 ' 

where, according to Eq. the a^ 1 ^ are given by 

= 9ii^ {£l£2} + 9i2al { ^\ (a = U,L) 
4^> = g^'^g^/^. (21) 

As before the curly bracket notation {£±£2} stands for either the unpolarized or the longi- 
tudinal spin-spin correlation case. Because a a ^ ili2 ^ and H a ^ xl ^ are proportional to each 
other, one can freely substitute H 3 ^ 1 ^ for cr^ 1 ^ in Eq. (|20|) . 

In the forward and backward direction (cos# = ±1) one finds the Born term relation 

P u (Born; cos 9 = ±1) = -1 (22) 

regardless of the energy y/q 1 which follows again from angular momentum conservation 
since the spins of the quark and antiquark have to be aligned in the forward and backward 
direction. Formally this can be seen by substituting the Born term relations f|T9l) into 
Eq. P}. 

Next we discuss the limiting values of the longitudinal spin-spin correlation asymmetry 
at threshold v — > and in the high energy limit v — > 1. At threshold v = one find|3 

threshold: P ee (Born; cos 9) = - cos 2 9 (23) 

which implies that the longitudinal spin-spin asymmetry vanishes at threshold at 90°, i.e. 
P u (Born; cos9 = 0) t h rC shoid = 0. The angular average is (P u (Born; cos 9) threshold) = — 1/3 
in agreement with Ref. [6]. In the high energy limit m q / 'y/q 1 — > where the quark and 
antiquark spins become completely aligned, i.e. o~l = crj^ l£2 ' > = 0, one has 

high energy limit: P ee (Born; cos 9) = —1 (24) 

independent of the polar orientation. 



2 Note that in the threshold region QCD binding effects modify the naive threshold results significantly. 
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4 0(a s ) tree graph contributions 



The 0(a s ) spin-dependent hadronic three-body tensor H fJil/ (pi,p2,P3, s[, s l 2 ) can easily be 
calculated from the relevant Feynman diagrams. In order to calculate the longitudinal spin- 
spin correlations one needs an explicit representation of the longitudinal spin components 
of the quark and antiquark. They are given by 



is 



ir 



^(V(i-y) 2 -£;0,o,i-y), (25) 



_ _f2 



£;(l-z)sin0i 2 ,O,(l-z)cos0i 2 ) (26) 

where 9\ 2 is the polar angle between the quark and the antiquarkU The polar angle is 
given by (£ = Am 2 Jq 2 = 1 — v 2 , y = 1 — 2p 1 g/ g 2 , z = 1 — 2p 2 g/ g 2 ) 



cos o a = , y z +y + z - 1+ t . (27) 



In the two-body Born term case, where y = z = 0, one has cos #12 = — 1. The polarization 
four-vectors (1251) and (|26|) then simplify to their expected two-body representations. 

For the spin-spin dependent piece one obtains (with N := a s N c CFq 2 /4:Tiv, R y : = 



y) 2 -£and R z :=J{\-zf-i) 



N 



R 3 yR Z 



2(2 -0(24- 23£ + 
1 



+2(8 - 90(1 - 0(2 - 0- + 4(30 - 24£ + 4£ 2 )y 

y 

-2(44 - 26£ + 3£ 2 )y 2 + 8(5 - £)y 3 - % 4 

+202 - 0(i - 2 (\ + \)- e^ 1 " 0(2 - 0-2 

\y z z z J z z 
2 34 

+402 - 2 ^ - 403 - 20^ + 202 - 0^ 

+2(16 - 110(1 - 0(2 - 0- " 4(1 - 2 (2 - 2 - 

z yz 

2 

-2(2 - 0(26 - 27£ + 30)- + 8(11 - 13f + 4£ 2 ) — 

z z 
3 4 

-4(1 -0(10-0- + 8(1-0- 

z z 



3 In the three-body case the quark and antiquark are no longer back-to-back in general and thus the 
spins in the aligned spin configurations (||) and (|t) are no longer parallel in general. However, in 
the (ii)-case the top and antitop quarks are quite "stiff" w.r.t. gluon radiation for the relevant lower 
energies. The average opening angle is (cos#i 2 ) = —0.996 and (cos#i 2 ) = —0.980 at yjcf = 500 GeV and 
^fq 1 = 1000 GeV, respectively [28]. 
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+8(5 - 30* -2£(1- 0(2-0 



-2(2 -0(6 -9, + e)- y -4(6 + 



+8(1 - 20^ - 4(2 - 2£ + £V + 2^(1 -0(2-0^ 
+4(2 - 0(1 - 20- + 8yz 2 - 16^yz 2 - 8y 2 z 2 

y 



H, 



U 



{tree) 



R 3 yR Z 



40 - 36£ - 2£(1 - 2 ^ - 2(8 - 90(1 - 0" 

1/ y 



-2(12 - Oy - 2(4 - i)y 2 + 8y 3 - 2£(1 - OS 

z z 

+6^(1 -04- 2^(3 -0^ + ^¥ 
-2(16 - 110(1 - 0- + 4(1 - 2 (2 - 0- 

2 

+2(24 - 23£ + i 2 )- - 2(16 - 50 — 

z z 

3 

+2(4 + 0- - 2(16 - 30* + 2^(1 - 04 
+2(8 - 11^ + + 2(8-0^ 

y 

+2(8 - O* 2 - 2^(1 - O^ " 4 (2 " 30- " 8^ 

3T 2/ 



H\^\tree) = 



N 



R 3 yR Z 



16 + 56e-43e 2 + 5e 3 -e 2 (l-0 2 A 

y 



-^(S - 90(1 - 0- + (32 - 70£ + 23e)y 

-2(8 - 90^ 2 - e^ 3 + 2& - ea - o 2 - 



2„,4 



+30(1 - 04 - 2£ 2 (2 - 0^ + H- - ^ 

z 2 z l 

-Oi6 - iiO(i - 0- + 20i - 2 (2 - 0- 

z 



yz 



+030 - 33^ + 50)- - 902 - 0- 
+02 - j - ^y- + (32 - 66£ + 23f 2 )* 



+O(i-o4+^( io - i5 e+3a- 



(48- 38^ + 0)^ 
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+2(8 - iWz + ^y 6 z - (16 - 26£ + OK - T (1 - 



-2£(3 - 40- + 2(8 + £)yz* + 2&/V 



H 



{tree) 



R 3 yR Z 



24 - 25£ + 30 - f (1 - 1 



1 



,J_ 



[8 - 90(1 - 0- " (22 - 90y + 2y 2 + 6y 3 - 2y 4 

y 



2Q/ 2 ty z | ty A 



z^ 



z^ 



-(16 - 110(1 - 0- + 2(1 - 2 (2 - — 

z yz 

+ (22 - 21£ + 0)- - (10 - 30— - (2 - 0— + — 

z z z z 

-13(2 - + oi - o4 + ( 10 - + 3 ^ 2 ) £ 

+ (18 + i)yz + 2y 2 2 - Ay 3 z + (10 + O^ 2 



^ 2 , . ^ 2 



2(3-40--2^ 2 -2yV 
2/ 



4(£i€ 2 ) 



(tree) 



2N 

rJr z 



32 + 26£ + e + 201 - 2 4 + 16 (! " O 2 - 



y 



+4(7 - £)y - 4(4 - Ol/ 2 + V + 201 - 0' 



(30) 



(31) 



-40l-Ojr + £(4-30^-02-0|r 

+6(4 - 30(1 - 0- - 4(1 - 2 (2 - 0- - 4^ 2 

z yz 

2 3 

-(28 - 32£ + 3£ 2 )| + 2(1 - 0(8 - 0^ - 4(1 - j 
+2(4 - 0(1 + 0* - 201 - 04 - (12 - 16£ + e 2A " 



+4(i - t)y Z - 4^ 2 + 20i - 0^ + (4 - - e)- 



(32) 



The unpolarized helicity structure functions are needed for normalization purposes. They 
read (see e.g. Refs. [H]) 



Hb(tree) = ^ 



„ a, ,,2 4 £ 60 20 16 32£ 120 
16 - 4^ - 4^ 2 - -| + -V - -V + 

y y y y y 



z A z A 



y 

2£ 3 { 8Q/ Aj 2 y Ajy 2 48 | 56£ 
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Hu(tree) 



Hl(tree) 



16£ 2 16 32£ 20£ 2 4£ 3 56y 28£y 2£ 2 y 

z yz yz yz yz z z z 
32y 2 8y 3 _ 8z 4£z 2£ 2 z 



+ — 8£z H 2 2 h 8?/^ 

z y y y 



, 2£ 2f 2 

+ 2^ - 4 + 4- 

y y 



o -|- Z£ - 



z c 



24 16f 

2; z 



12£_ 2^ + 2^ 

y z 2 z 2 

r2 



2/2 

2£z 



+ 

yz yz 



24y 2£y 8y 2 
H — — + 8z + 

z z z y 



N 



Hl(tree) = — 



y y y y 



R 2 

_ J e + e + zey + ey 



IQy + 4£y + 4£y 2 



2„,2 



| 2jy | 2£ 2 y | 2{jy* 



£2 c3 
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(33) 



(34) 
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(37) 



5 One-loop contributions and 0(a s ) corrections 



The one-loop contributions to the unpolarized helicity structure functions can be taken 
from Refs. [T3l [14] . Let us first list the contributions from the real part of the one- loop 
amplitude. For the unpolarized case one has 



Hhiloop) 
H 2 (loop) 
Hl(loop) 
H 4 F {loop) 



AN c q 2 (ReA + v 2 ReC), 

AN c q 2 (ReA-v 2 ReC), 

Hl(loop) = 2N c q 2 (£ReA + v 2 ReB) 1 

4N c q 2 v(ReA + ReC). 
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For the longitudinal spin-spin correlation components one obtains 
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L 
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-AN c q 2 v(ReA + ReC). 



(39) 



Since the one-loop contributions refer to two-body final states they also satisfy the two- 
body relations ( TT9l . 

There are also contributions coming from the imaginary part of the vertex correction 
which multiplies the imaginary part of the Breit-Wigner function of the Z resonance as 
indicated in the electroweak model parameters and gi3 in Eqs. ([7]) and ffl^l) . The relevant 
hadron tensor component results from the V/A interference term in the F projection. 
Although these contributions are rather small (especially when one is far away from the Z 
resonance) they are included in our numerical results. The imaginary part contributions 
are given by 
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Hl(loop) 
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The one-loop form factors A, B, and C appearing in Eqs. fl38l) and fl39|) read (Cp 

4/3) [n 
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(41) 



The one-loop contributions have been infrared regularized by introducing a (small) gluon 
mass m g = y/Aq 1 following Refs. [H] . 

What remains to be done is to perform the phase space integrations over the tree-graph 
contributions listed in Sec. 4. One first integrates over z and then over y. The relevant 
integration limits incorporating the auxiliary gluon mass are 
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(42) 



and 



A£ + A, y + = l-J£. 



(43) 



13 



The introduction of a (small) gluon mass distorts phase space away from the infrared 
singularity at y = z = 0. The infrared singularities in the tree graph and one- loop con- 
tributions cancel and one remains with finite remainders. It is quite clear that the finite 
result in the sum is independent of the specific regularization procedure. Finally, adding 
in the above one-loop contributions one obtains the 0(a s ) corrections (loop + tree). For 
the sake of completeness we include in our results also the unpolarized hadron tensor com- 
ponents which are needed for the normalization of the longitudinal spin-spin asymmetry. 
One obtains 
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(54) 



The unpolarized hadron tensor components Hjj(a s ), H^(a s ), H\{a s ), H^(a s ), and Hp(a s ) 
including the 0(a s ) integrated rate functions U (i = 1,2,3,6,7,8,9,10,12) have been 
calculated and listed before in Ref. [H]. The integrated rate functions £, (i = 13, 14, 15, 16) 
appear in the evaluation of the longitudinal spin-spin component H^ 1 ^ and have been 
listed in Ref. [6]. In addition to the decay rate functions calculated in Ref. [HI [33] the 
spin-spin polar dependence contributions bring in a set of new decay rate functions £j 
(i = 17, 18, 19, 20, 21, 22). The additional set of decay rate functions needed in the present 
application is given by 
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where 
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(61) 



We have checked that the sum of H]f 2 \a s ) and iJ^' 2 ^ l£2 ^(« s ) agrees with the form 
Hu+L, li2 \ a s) calculated in [Bj. 

6 Massless QCD and the zero-mass limit of QCD 

Before we turn to the numerical evaluation of the longitudinal spin-spin correlations, we 
would like to discuss the m q —>■ limit of our analytical results with the aim to compare 
them to the corresponding results obtained in massless QCD (m q = 0). As is well-known 
by now, the m q — > limit in particular of the spin-flip contribution does not coincide with 
that of m q = QCD where there is no spin-flip US UHl HH H31 HI]- 

The m q = expressions can be calculated in dimensional regularization as described in 
Ref. [14]. For the 0(a s ) contributions one obtains 

rrl/ \ a at 2 Q sCf 1(^ 2 ), . 2 a sC-F 

H u( a s) = 4N c q , H v y \a s ) = -AN c q , 
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The 0(a s ) m q = single-spin functions ifa^ 1,2 ^(a; s ) are given in Ref. [H] or can be read 
off from Eq. ( 16S1 . 

Let us now consider the m q — > limit of the integrated decay rate functions ti discussed 
in Sec. 5. One obtains 
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Using these limiting expressions one obtains for the 0(a s ) m q 
functions (loop + tree) 



(63) 

unpolarized structure 
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and, for the longitudinal spin-spin correlation functions, 
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With the square bracket notation we indicate the difference between the 0(a s ) m q — > and 
m q = results which we call anomalous terms@ Note that the current-current structures 
for j — 2,3 are zero in the limit m q —>■ (and also for m q = 0). 



4 We use the phrase "anomalous flip contribution" since the same anomalous flip contribution contributes 
to the absorptive part of the VVA triangle diagram in the zero mass limit of QCD which in turn can be 
related to the well-known axial anomaly via a dispersion relation approach [29j . 
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For completeness we also list the 0(a s ) m q — > single-spin functions using the same 
square bracket notation. They read [H] 
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We finally collect all our m q — > results, including the Born term contributions 
Eqs. (jl~6j) and (ITT)) , where we make use of the representation Eq. ([3]). Again we split 
off the anomalous contributions using the square bracket notation. One obtains (Cp = 4/3 
is made explicit here) 
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In Tables [JJ and [2] we list all m g = and m q — » contributions to the various spin 
configurations for the parity-even (VV) and parity-odd (VA) current contributions using 
the representation Eq. (ISTj) . We have again used the square bracket notation to denote the 
anomalous contributions. The nonvanishing anomalous spin-flip contribution proportional 
to [4/3] in Tabled] e.g. in the (||) spin configuration agrees with the corresponding result 
in [16] where this contribution was referred to as the chirality breaking contribution. 
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Table 1: Born term and 0(a s ) corrections to specific spin configurations (first column) in 
QCD(m g = 0) and QCD(m g — > 0) for the parity-even current contributions. The entries are 
given in terms of contributions to the hadron tensor components H^ v (s{ , s|) = H^ A (s{, s l 2 ) 
(a = U, L, F) in units of N c q 2 . Anomalous contributions are shown in square brackets. 



Using the m q — > limiting expressions Eqs. (164")) and (|65|) one finally obtains for the 
longitudinal spin-spin asymmetry 



a s C F (^ , / (1 + cos 2 6> + 4 sin 2 %n \ a s C f 
(cos 6^) = — 1 + |4J- 



4tt \ \(1 + cos 2 %n + 2cos9g 44 J 4tt / 
= -l + [4]^ + 0(a 2 ). (68) 

where we again have encased the "4" in square brackets in order to identify the contribution 
as anomalous. The result in Eq. (1681) is independent of the flavour of the quark and thereby 
also holds for polarized lepton pair production e + e~ — > fi + fi~,T + r~ replacing, of course, 
C F a s by a. We emphasize again that the normal m q = contributions lead to the value 
P (cos#) = —1 in the high energy limit. It is quite remarkable that the anomalous 
contributions do not show any cos 6* dependence at order 0(a s ). 



Anomalous contributions from the 
universal splitting function 



It was shown e.g. in Ref. [18] that, in the limit m e — > 0, there is a finite helicity flip 
radiation (i.e. — > en + 7 and — > + 7) which arises from collinear photon emission. 
It can be described in terms of a universal splitting function -D|j om = f-z where z is the 
fractional energy of the photon (not to be confused with the energy-type variable used in 
this paper). In this order of perturbation theory the derivation of Ref. [18] can be directly 
transscribed to the QCD case where one now has 



D ™ = C ^ z - (69) 
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Table 2: Born term and 0(a s ) corrections to specific spin configurations (first column) 
in QCD(m g = 0) and QCD(m g — > 0) for the parity-odd current contributions. The en- 
tries are given in terms of contributions to the hadron tensor components Hy A (s{, s^) = 
H AV {s\, si), Hl A (s{, s e 2 ) = H AV (s{, si), and H^ A (s{, si) = H AV (s{, si) in units of N c q\ 
Anomalous contributions are shown in square brackets. 



Integrating over the fractional energy z gives 

rl ry ry 

H~= / H{Born)C F —zdz = C F -^H(Born) (70) 
1 Jo 2ir An 

where we have, for the sake of simplicity, dropped all indices on the hadron tensor H(Born). 

One knows that there is no anomalous contribution to a given fermion transition when 
one integrates over the whole two-body phase space and one sums over the two spin states 
of a given final fermion. This was explicitly demonstrated for muon decay — ► e~ u e 
in Ref. [19] where the 0(a) contributions to the two final electron spin states — > e# + 7 
and cl —>■ cl + 7 were calculated in the limit m e — > 0. There was an anomalous spin-flip 
contribution to cl — ^ &r + 7 determined by the universal splitting function in Eq. (f59"j) 
accompanied by an anomalous non-flip contribution toei^ei + 7of equal and opposite 
strength. When summing the spin-flip and non-flip contributions the anomalous terms 
cancelled leading to a "normal" unpolarized rate. This implies that, when integrating 
over the whole phase space, every anomalous spin-flip contribution is accompanied by 
an anomalous non-flip contribution of equal and opposite strength. As shown e.g. in 
Ref. (18] the photon in the collinear emission process e F — > en + 7 is lefthanded, i.e. 
one has e F — > &r + 7l in agreement with angular momentum conservation. By the same 
token the anomalous non-flip contribution must vanish in the forward direction as has been 
explicitly demonstrated in [18J. 

The same empirical pattern is observed in the m q —>■ limit of e + e~ — > qqg. Consider, 
for example, gluon emission from the quark and keep the antiquark helicity fixed at q^. 
The anomalous spin-flip contribution q F — > qn + g is determined by the universal splitting 
function Eq. (1691) which is cancelled by the residual non-flip contribution qi — > qi + g when 
one takes the sum of the two, i.e. the sum has no residual mass effects. We have verified by 



21 



explicit calculation that the residual non-flip contribution originates from the near-forward 
region irrespective of the fact that the non-flip contribution vanishes in the exact forward 
direction. 

In this way one can determine the anomalous flip and non-flip contributions to each 
spin configuration. In fact, one reproduces the anomalous entries in Table [Hand Table [2] 
by using the relations 
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and correspondingly for the spin configurations [||] and [jj]. The anomalous contributions 
to the spin-spin structure function Hy 2 can be seen to be 50% flip and 50% non-flip. The 
anomalous contribution to the single-spin structure function derives entirely from the 
anomalous non-flip contribution contrary to what was stated in [14] . This, in particular, 
means that the universal anomalous spin-flip contribution will not show up in single-spin 
polarization measurements, but contributes at the 50% level to spin-spin polarization mea- 
surements. 

One can explicitly verify from the entries of Table [Hand Table [2] that the anomalous flip 
and non-flip contributions to a given quark or antiquark transition cancel upon summation. 



We mention that the "normal" 0(a s 
determined from Eq. ([6 



contributions in Table [T] and Table [2] can be 



8 Numerical results 

We are now in the position to discuss our numerical results for the cos ^-dependent nor- 
malized longitudinal spin-spin asymmetry P (cos 9) defined in Eq. (1201) . 

In Fig. [T] we show the cos 9 dependence of the 0(a s ) longitudinal spin-spin correlation 
function P u (cos8) in e + e~ — > ti for the four different center-of-mass energies \fq* = 360, 
400, 500, and 1000 GeV. At y/q 1 = 360 GeV one is sufficiently far above the (tt) -threshold 
for perturbative QCD to apply. The above range of center-of-mass energies is the envis- 
aged range for the proposed linear colliders. At the lowest shown energy y/q* = 360 GeV 
the cos 9 distribution of P u is still quite close to the threshold Born term distribution 
P ee (cos9) = — cos 2 9 (see Eq. (1231)). However, one already notes a small forward-backward 
asymmetric effect which is induced by the (^/^-interference term Hp i2 ^. In the forward 
and backward directions cos 6* = ±1 the longitudinal spins are quite close to being 100% 
aligned as is true for the Born term case (see Eq. f[22|) ). Away from the forward and back- 
ward direction the destructive longitudinal contribution H^ 1 ^ (j = 1, 2) comes into play 
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Figure 1: 0(a s ) polar angle dependence of the longitudinal spin-spin correlation asymmetry 
in e+e" -> it at V^F = 360, 400, 500, and 1000 GeV. We always take m t = 175 GeV as 
default value. 
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k tter. At ^q 1 = 1000 GeV 
75% over most of the range of 



creased, H^ 1112 ' becomes smaller and thus P u (cos9) becomes flatter. At ^Jq 1 = 1000 GeV 
the alignment of the spins of the top and antitop exceeds 
cos 9. 

In Fig. [2] we exhibit the relative size of the 0(a s ) corrections which we define by the 
measure 
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Fig. [2] shows that the 0(a s ) corrections to P ee (cos8) are strongly polar angle dependent. 
They are not very small and can reach values as high as ~ 4%. It is quite remarkable that 
the radiative corrections for the lowest energy y/q 1 = 360 GeV are very close to zero in 
the forward and backward directions. As the energy is increasing the radiative corrections 
in the forward and backward directions become larger and they start approaching their 
anomalous limiting value Eq. fl68|) . 

In Fig. [2] we have also included a plot of AP u (cos8) for the rather large center-of-mass 
energy of yfq 1 = 10 4 GeV (we are not implying that such center-of-mass energies will ever be 
reached in terrestrial laboratories). The purpose is to exhibit the anomalous contributions 
induced by the 0(a s ) corrections when m q /y/q^ —* 0. In terms of the measure AP u (cos8) 
defined in Eq. (173|) one has in the high energy limit (see Eq. (l68|) ) 
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Figure 2: Relative size of 0(a s ) corrections to the longitudinal spin-spin correlation asym- 
metry in e+e" -> tt at y/q* = 360, 500, 1000, and 10000 GeV 
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Figure 3: Relative size of 0(a s ) corrections to the longitudinal spin-spin correlation asym- 
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Using a s {\f(p = 10 4 GeV) = 0.0712 one numerically has AP ee (cos6) = —0.0303 which is 
quite close to the value ~ —0.026 shown in Fig. [2]. The remaining difference can be traced 
to non-asymptotic effects. 

In order to highlight the importance of hard gluon emission for the radiatively cor- 
rected longitudinal spin-spin asymmetry we have determined AP u (cos9) in the soft gluon 
approximation as has also been employed in Ref. [3]. Fig. [3] shows a plot of AP ee (cos6) 
in the soft gluon approximation with the gluon energy integrated to its maximal energy. 
Since the soft gluon approximation is essentially a two-body approximation the radiative 
corrections can be seen to vanish in the forward and backward directions as discussed in 
Sec. 3. For the lowest energy yfq 1 = 360 GeV the corresponding two curves in Figs. [2] 
and [3] practically lie on top of each other since there there is very little phase space left 
for gluon emission. At yfq 1 = 500 GeV the difference between the full and approximate 
calculations becomes noticeable. In particular the cos ^-dependence of the full calculation 
is more pronounced. The same holds true for ^/q I = 1000 GeV where the difference be- 
tween the full and approximate calculation becomes as large as ~ 1% in the forward and 
backward directions. At \fq^ = 10000 GeV one has AP u (cos8) ~ in the the soft gluon 
approximation as expected whereas AP^(cos#) ~ —2.6% in the full calculation due to the 
anomalous contributions as discussed before. 

The relatively small value of the 0(a s ) correction for the longitudinal spin-spin asym- 
metry in Fig. [2] has to be contrasted with the large 0(a s ) corrections to the cross section as 
shown in Fig. HI For example, at y/q 1 = 500 GeV the 0(a s ) corrections to the cross section 
amount to 15.4%. The reason that the radiative corrections to P u are smaller than those 
for the rate is that the 0(a s ) corrections in the numerator and denominator of P tend to 
go in the same direction and thus tend to cancel when one takes the ratio. However, this 
is not true for the anomalous contribution which contributes only to the numerator. 

It is quite remarkable that the radiative corrections to P u in the forward and backward 
directions are not small for the higher energies even if they only arise from non-Born- 
term-like hard gluon emission. Compare this to the corresponding case of top quark decay 
t — > b + Wr where the right-handed Wr is again only populated by non-Born-term-like 
hard gluon emission. In the case the radiative corrections amount to only 0.1% of the total 
Born term rate J3U|. The explanation is that in the e + e~ — > ti case one is sensitive to the 
enhanced forward hard gluon emission region whereas in the t — > b + Wr case there is no 
such enhancement. 

In Fig. Owe display the polar angle dependence of the 0(a s ) differential cross section 
for (tt)-production in e + e~-annihilation. As Fig. shows, the cross section strongly peaks 
in the forward direction. For example, for yfq 1 = 1000 GeV there is more than a factor of 
10 cross section difference between the forward and backward directions. The fact that the 
cross sections peak so strongly in the forward direction is quite welcome from the point of 
view of checking the 0(a s ) prediction for the longitudinal spin-spin asymmetry which, as 
remarked on before, is close to its m q = value P u = —100% in the forward direction. 

In Figs. EHH] we show the corresponding plots for (66)-production at center-of-mass 
energies of \/q* = 10, 20 and 50 GeV and at the Z pole (mz = 91.188 GeV). We use a pole 
mass of mb = 4.8 GeV for the bottom quark. In Fig. [6] we show the polar angle dependence 
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of the 0(a s ) longitudinal spin-spin correlation function P u (cos6). The forward-backward 
asymmetry of P u is not as pronounced as in the top quark case and, contrary to the top 
quark case, its angular distribution is only slightly enhanced in the forward hemisphere. At 
the Z pole one has an almost flat distribution close to the value of P = — 1 as expected 
for massless QCD. The deviation from the —100% value is again due to the anomalous 
contribution. 

In Fig. [7] we exhibit the size of the 0(a s ) corrections defined by the measure AP U as 
defined in Eq. ( 1731) . Close to threshold at y/q 1 = 10 GeV the 0(a s ) corrections in the 
forward and the backward directions are quite small. The corrections become larger away 
from the forward and backward directions and rise to ~ 3.5% at 9 = 90°. As the energy 
increases, the 0(a s ) corrections become progressively flatter and tend to turn to negative 
values. At the Z peak the 0(a s ) corrections are almost flat and reach a value of ~ —3.4% 
which again shows that they are dominated by the anomalous contribution. 

In Fig. [S] we show the cos 6*-dependence of the (66)-production cross section. The distri- 
butions are much flatter than in the top quark pair production case. The differential cross 
section falls as the center-of-mass energy increases and then reaches back to high values 
on the Z peak. 



9 Summary and conclusion 

We have presented analytical results for the 0(a s ) corrections to the polar angle dependent 
longitudinal spin-spin asymmetry P ee (cos6) in e + e~ — > qq. A numerical evaluation of 
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P u (cos6) for top quark pair production shows that the 0(a s ) corrections are not small and, 
depending on the energy, can be strongly polar angle dependent. The 0(a s ) corrections 
to P are smaller than the 0(a s ) corrections to the cross section itself. The reason is 
that the 0(a s ) corrections to the numerator and the denominator in the (normalized) 
longitudinal spin-spin asymmetry tend to go in the same direction with the effect that the 
0(a s ) correction to the ratio becomes reduced. At lower energies the longitudinal spins 
of the quark and antiquark become almost 100% aligned in the forward and backward 
direction, i.e. one has P u (cos8 = ±1) ~ —1 with small radiative corrections. As the 
energy is increased the radiative corrections become larger at the endpoints signalling a 
slow approach to the asymptotic regime. 

For bottom quark pair production on the Z peak one is already close to the asymptotic 
regime since m^/m| = 2.8 x 10~ 3 (using nif, = 4.8 GeV, rriz = 91.188 GeV as before). In 
massless QCD one expects P^(cos 8) = — 1 to any order of a s because there are no spin-flip 
contributions in massless QCD. This is different in the m q — > (or high energy) limit of 
QCD where there are residual mass effects starting at 0(a s ) which change the naive pattern 
of massless QCD. In fact, in the high energy limit (or for m q — > 0) one obtains instead 
P ee (cos8) = — 1 + 4a s /3ir due to the anomalous 0(a s ) flip and non-flip contributions. 
The numerical evaluation of the longitudinal spin-spin asymmetry for bottom quark pair 
production on the Z peak shows that one is already close to the asymptotic form but 
noticeable preasymptotic effects are still present with respect to the flatness of the cos 8 
distribution and the size of the anomalous contributions. We have discussed in some detail 
the role of the anomalous contributions and how they contribute to spin-flip and non-flip 
terms in the three correlation structure functions that determine the cos ^-dependence of 
the longitudinal spin-spin asymmetry. 

All our results can be applied to QED final state corrections in polarized lepton pair 
production in e + e~ — > fi + fi~ and e + e~ — > t + t~. The radiative corrections are, of course, 
smaller in the QED case due to the smallness of the QED coupling constant a. In particular, 
in the asymptotic limit, where P u (cosd) = — 1 + a/ir, the anomalous contribution to the 
longitudinal spin-spin asymmetry is now only of 0(0.23%). 

What remains to be done is to calculate the radiative corrections to the remaining 
transverse-longitudinal and transverse-transverse spin-spin density matrix elements in the 
helicity system. Alternatively, one can study the spin-spin density matrix in the off- 
diagonal basis representation of Ref. [2], which is optimal at the Born term level. Radiative 
corrections to the spin-spin density matrix in the off-diagonal basis have been considered 
before using the soft gluon approximation [3J. It is, however, clear that, in the soft gluon 
approximation, one misses potentially large contributions from e.g. the near-forward emis- 
sion of hard gluons as described in this paper. 
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